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The properties of the s-wave for a quasi-free particle with position-dependent mass(PDM) 
have been discussed in details. Differed from the system with constant mass in which the 
localization of the s-wave for the free quantum particle around the origin only occurs in two 



in 
o 
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dimensions, the quasi-free particle with PDM can experience attractive forces in D dimen 



sions except D — 1 when its mass function satisfies some conditions. The effective mass of 



, a particle varying with its position can induce effective interaction which may be attractive 

in some cases. The analytical expressions of the eigenfunctions and the corresponding prob- 
ability densities for the s-waves of the two- and three-dimensional systems with a special 



' PDM are given, and the existences of localization around the origin for these systems are 

. shown. 

o . 

^ . PACS numbers: 03.65. Ca, 03.65.Ge,02.30.Hq 

o 

\Q ■ Keywords: position-dependent mass, s-wave, Schrodinger equation, dimensionality of the space, 



localization, eigenfunction, probability density 



INTRODUCTION 



The properties of a system interested in physics have a close relationship with its dimension- 
ality of space. Quantum Hall effect is a macroscopical quantum phenomenon occurred in the 
two-dimensional system Bose-Einstein condensation only takes place in the three-dimensional 
Bose systems but not in one- and two-dimensional casesP]. There is no finite-temperature phase 
transition in the one-dimensional system, and there may exist Kosterlitz-Thouless transition of the 
topological type in two-dimensions jjj and the conventional thermodynamics phase transition in 
three-dimensions Q|. One has made extensive theoretical and experimental researches on properties 
of the systems with various dimensions and also developed many effective methods to understand 
them. So far, quantum mechanics is believed to be the most successful theory to tackle various 
problems in physics. 

In quantum mechanics, we usually solve Schrodinger equations of various potentials to obtain 
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the eigenfunctions and the related physical quantities of the system. It is well known that the wave 
functions for both the two-dimensional system with axial symmetries and the three-dimensional 
system with spherical symmetry can be separated into the product of the radial and angular parts. 
Due to the fact that the angular component of the eigenfunction is the spherical harmonics(in three- 
dimensions) or its reduced form(in two-dimensions), the emphasis to solve Schrodinger equation 
is put on its radial component for the system with spherical or axial symmetries 4] . The radial 
Schrodinger equations for two- and three-dimensional systems have little difference in their forms, 
we can transform one into another by just making a simple replacement of the quantum number 
of angular momentum. It is in this sense that there is no remarkable difference between two- and 
three-dimensional systems in the framework of quantum mechanics as in other problems. However, 
Cirone and coworkers have made systematic studies on the two-dimensional system through solving 
the Schrodinger equation in recent years P. Ifl I7 . P. Ifll llO . 11, 3|. They find that the s-state of a 
free quantum particle in two dimensions can have a unique localization around the origin which 
does not present in other dimensions. 

The study of quantum systems with position-dependent mass(PDM) has become one of the 
active subjects over the last years. These systems have found wide applications in the determination 
of electrical properties of various systems such as semiconductor, quantum dot, liquid crystal and 



so on. Theoretically, the search for the exact solutions of the wave equations, especially t 



e wave equations, especially those of the 

Ea, 14, 13, C 12, 3, Q, , 21, 22, 3 



Schrod ing er equation with PDM has aroused great interest! 
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, they may provide useful models to the description of the systems with PDM. 
A natural problem is what relations there exist between the dimensionality of space and properties 
of the system in the case of PDM. Up to now, the study of the system with PDM is mainly focused 
on one rather than higher dimensions. As far as we know, there are no investigations on the possible 
differences between systems of different dimensions with PDM. In this article, we will discuss the 
properties of systems of different dimensions with PDM through their zero angular momentum 
states (s-state). We find that the effective attractive potential appeared in the two-dimensional 
system with constant mass also exists in the two-dimensional system with PDM, and in particular 
that there is an additional effective potential related to the mass function of the particle in the 
latter case. The effective potential induced by the variation of mass with position may be attractive 
or repulsive depending on the sign of the first order derivative of the mass function with respect 
to the position. The more important thing is that this induced effective potential appears in the 
system with any dimension other than one. When effective mass is a constant, the induced effective 



potential vanishes, which reflects the fact that the system with PDM is more general than that of 
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the constant mass. In this article, we will give the explicit expression of the effective potential for 
the quasi- free particle with PDM in D dimensions. For the so-called quasi- free particle we mean 
that the particle does not subject to other potential field except its mass is position-dependent. 
With the mass being constant, the quasi-free particle becomes free one. We also solve the radial 
Schrodinger equation to get the eigenfunctions of the s states and the corresponding probability 
density functions, and analyze their behaviors for a special mass function in the two- and three- 
dimensional systems, respectively. 

This article is organized as follows: In section [HI we will analyze the s-state of the quasi- free 
particle with PDM and make a comparison between the properties of the systems in different 
dimensionality D as well as btween those of the system with constant mass and with PDM in 
terms of the effective potential. In sections IIII Al and IIIIBI we will get the exact solutions for 
the s-wave of the quasi-free particle with a special mass function in two and three dimensions, 
respectively. They localize around the origin, which show the existence of the effective attractive 
potential. In Section llVl we will make some conclusions and discussions. 

II. RADIAL SCHRODINGER EQUATIONS AND THE EFFECTIVE POTENTIALS OF 
D-DIMENSIONAL QUASI-FREE PARTICLE WITH PDM 

We first consider a D-dimensional system with PDM. We will assume that the effective mass 
varies with the radial coordinate m = m(r). When the mass of a particle depends on its position, 
the mass and momentum operators no longer commute, so there are several ways to define the 
kinetic energy operator of the quantum system. In this paper, we will adopt the form of the kinetic 
energy introduced by Levy-Leblond and etcQ], which reads 

T-— p— ^- p (1) 

where /x(r) = is dimensionless mass, and the momentum operator p can be expressed in terms 
of L>-dimensional differential operator Vb as 

p = —iKS/d- 

In Cartesian coordinate, we have the representation of the operator V_d to be of the form Vz> = 
(^ii ' ' ' ' ' ' dxD^' ^ n ^ e following sections, we will use the units of mo = l,h = 1, then the 
Schrodinger equation of the quasi-free particle with PDM in D-dimensional space may be written 



4 



-V — 

2 D W) 



Eip(r). 



We adopt spherical coordinate (r, Oi, 62, ■ ■ ■ , Od-2, 4>)i an d define a operator 
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£>r, 
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where 



1 



(r ■ p + p ■ r) 



1 



.(D-l)/2 



is the radial momentum operator in D-dimensional space 



d D-l 
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The operator 



T r reduces to t£ with the mass being constant. We now rewrite the Hamiltonian of the quasi-free 



particle with PDM in D dimensions as follows 

1 
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H 



Mr) 



L 2 , I (D-l) (£-!)(£> -3) 



2r 



4r 2 



(4) 



where L is the D-dimensional angular momentum operator. 

For the ZMi— , SyS te m with .WflUB^ the ei — io„ , W of the SyS ,e m 
can be made separation of variables 



ran spnerical symi 
V;(f) = r-^-^^Rir) • 



/ 1 l f7 l; t7 2, 



,0D- 



(5) 



where 2 ... ^(#1, #2, • • ■ ,0d-2i<P) is D-dimensional hyperspherical harmonic function, and / is 
the quantum number of angular momentum. 

With Eq.Q and the eigenfunction we can get the D-dimensional radial Schrbdinger equa- 
tion from Eq.@ 

d 1 



dr 2/j,(r) dr 2[i(r 

1 ia+p-2) 



d 1 1(1 + D -2) // (D-l) (D-l)(D-3) 



4r 



8/ir 2 



i?(r) = ER(r) 



(6) 



where 2//( r ) H '~^ Z-> * s ^ ne centrifugal potential which results from the action of the term jj^r)!^ ™ 
Eq.@ on the hyperspherical harmonic function in Eq.(JSJ). The factor 1(1 + D — 2) is the eigenvalue 



of the operator L 2 with eigenfunction being the hyperspherical harmonic function 



I . If we 



only consider s-wave, the contribution from the centrifugal potential is zero. The last two terms 
in the bracket on LHS of Eq.@ are the quantum effective potential (QEP) for the system with 
PDM, 



Vqep 



VHP-1) , (D-l)(D-3) 



4r 



+ 



(7) 
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The first term — Ay ^r"^ m Eq.0 is the effective potential due to the dependence of the mass on 
position;The second term comes from the reduction of the motion of the particle in D dimensions 
to that of radial component and it also appears in the system with constant mass. 

Now we can make a few remarks on Eq. (JTJ) : (i) A quasi- free particle with PDM is generally not 
free in the radial direction. When D = 3 and the mass of the particle is constant, QEP disappears. 
However, in the system with PDM, the potential —2^27 originated from the variation of mass with 
position does not vanish. For D = 1, QEP is zero for systems with constant mass or PDM. So 
QEP remains to be related to the dimensionality of the system, (ii) For D > 3, QEP is a repulsive 
potential if the mass of the particle decreases with the radial coordinate, i.e. //(r) < 0. (iii) Being 
different from the system with constant mass in which attractive QEP only occurs for D = 2, we 
may also have attractive QEP for other D in the systems with PDM. The sign of QEP in Eq. 
determines whether QEP is attractive or repulsive. For example, the QEP is attractive for D = 3 
when the mass function p{r) satisfies the condition fi'(r) > 0. 



III. EIGENFUNCTIONS AND LOCALIZATION OF THE s STATES FOR A 
QUASI-FREE PARTICLE WITH PDM 

A. Two-dimensional system 

When D = 2, the Hamiltonian of a quasi-free particle with PDM in polar coordinate (p, (p) can 
be written as follows from Eq.Q 

where p p = —ifr{-§2 + ^) is the radial momentum operator for the two-dimensional system. 

The wave function of the two-dimensional system with axial symmetry can have the following 
form 

m = \r{pV 1ip , (9) 

VP 

where I = 0, ±1, ±2, • • • is the quantum number of angular momentum. If we only consider s-state, 
that is / = in Eq.©, then the second term in Eq.Q will disappear in the radial Schrodinger 
equation. The last two terms in Eq. (jHJ) is the QEP mentioned-above and is the reduction of Eq.Q 
in two dimensions, 

Vqfp = "M^F " WWp (10) 
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From Ea.(|10|). the following remarks are in order: (i) When mass is a constant, i.e. p! = 0, Vqep 
reduced to that studied by Cirone and coworkers (ii) If the mass p 

increases with r, that is fjf > 0, then Vqfp is always negative, and so it is a central attractive 
potential. In this case, the variation of mass strengthens the attractive force, (hi) If p! < 0, QEP 
may be attractive or repulsive depending on the specific form of the mass function. However, the 
dependence of mass on the position may weaken and even destroy the attractive action in this 
case. So in the two-dimensional system, a quasi-free particle with PDM is subject to an effective 
potential in its radial motion. This effective potential has an additional part from the variation of 
effective mass except for that occurred in the system with constant mass. 

Now we devote ourself to the solution of the s-state and discuss its properties for a given mass 
function. We have known that if the effective mass p increases with r, i.e. p! > 0, then QEP 
Vqfp is always attractive. When we only consider s-wave, the effect from the centrifugal potential 
is suppressed, so we expect that there exists the bound state localized around the origin and the 
corresponding probability function displays a maximum near the origin. 

We chose the mass function to be of the form p(p) = p a , and we also assume that the energy 
be negative, i.e. E = — \E\ so that the state is bound. Notice that the mass function may have a 
positive multiple factor, we group it into m®. The requirement that the mass function satisfies the 
condition pi > imposes the restriction a > on its exponent a. This is a stronger restriction for 
a. Because we only require that Vqfp be attractive, that is Vqfp < 0, Eq. implies that a > — ^ 
is sufficient for that requirement with the above mass function. Due to the following considerations, 
we need not worry about the possibility of the above effective mass increasing infinitely: We will 
soon see that the eigenfunction decays exponentially with the distance from the origin and so the 
increment of the effective mass is finite around the origin. On the other hand, we mainly concern 
the behavior of the eigenfunction for the s state around the origin. This situation is similar to that 
of a particle under the harmonic oscillator potential. 

Now the radial Schrodinger equation for the quasi-free particle with effective mass of the form 
p(p) = p a can be written as 

1 d 2 u(p) 1 — a du(p) I 2 , , , v 

where u(p) is the radial wave function which is related to R(p) in Eq.® through the equation 
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For s-wave, i.e. I = 0, Eq. reduces to the following equation 

d 2 u 1 — a du(p) 



dp 2 p dp 

Setting 



2\E\p a u(p) = 0. (13) 



z = 2^ (ot+ . a)/a y = a = 
a + 2 2 + a 



we may transform Eq. (|13|) into the standard modified Bessel equation 



36] 



^ + ^-(^ + ^W-a (15) 

Ea.p5|) has two special solutions I v {z) and K v {z) which are the first and second modified Bessel 
function, respectively. It should be noted that the range of is — | < ^ < 1 due to a > — |. 

We require that u be finite when p is large so the eigenfunction may be normalized. This is 
equivalent to £(z) — * when z — ► oo. Considering of the asymptotic behavior of the functions 



36j, we have the general solution of 



l v and K v with the modulus of their arguments being large 
Ea.(|15[) which is finite with large |z|, 

i = cK v {z\ (16) 
where c is a arbitrary constant. So the solution of Eq.(|13|) is 



u(p) = cz v K u (z) - ... 



^[2\E\ { l-v)]^z»K u (z) 



2\E\ sin(i / 7r)(l — v\ ,, . r— — i 



(17) 



UTT 2 



We can show that this eigenfunction has been normalized 

+oo r+oo r2ir r+oo 

dx dy\ip(x,y)\ 2 = dip pdp\u(p)\' 

oo J — oo J J 



2 sm(yn) 



roo 

\ zK 2 (z)dz = l. 
Jo 



(18) 



The probability density that the particle may appear between p and p + dp is 

f-2TT 

W(p)dp = / d9pdp\ip\ 2 = 2K\u{p)\ 2 pdp 



o 



4I.EI sin(^7r)(l — u) it*.,,,., , / . =n , 



(19) 



Z/7T 

In Fig^ we depict the probability density function W(p) with = i and ^ = —0.2,0.1,0.4,0.7, 
respectively. These figures indicate that W(p) sharply collapses at the origin and decays exponen- 
tially for large p, with a maximum close to the origin. In addition, the maximum of W deviates 
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from the origin as increasing v. In order to show how W(p) varies with p, we give the cuts of 
FigE along the radial in FigEl It is seen that the variation of W(p) with p is similar to that with 
constant mass as given in For other energies E, the corresponding W(p) has similar behavior, 
they all have the localization around the origin. 

The localization of the probability density for the s state of a quasi-free particle with PDM is 
produced by above QEP Vqfp. Because the mass function is chosen so that p! is always positive, 
Vqfp is an attractive potential. The above localization of the probability density confirms that 
the potential is attractive. 




(c) v = 0.4 (d) v = 0.7 

FIG. 1: The radial probability density of s-state of a two-dimensional quasi-free particle with PDM is at 
the origin, but has a maximum close to the origin. The radial probability density decays exponentially for 
large p. The parameters are \E\ = \, and v = —0.2, 0.1, 0.4, 0.7, respectively. 
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FIG. 2: The cut of the radial probability density of the s-state for a two-dimensional quasi-free particle with 
PDM shows that the maximum of the probability density deviates from he origin with v increasing. The 
parameters are the same as those in Fig^ 

If E = \E\ > 0, we make the similar transformation as in Eq. (|14l) . then the radial Schrodinger 
equation for / = and the mass function p{p) = p a reads 

az z dz 

Eq. (|2Uj) has two special solutions J v {z) and Y u (z) which are first and second kind of Bessel functions, 



respectively 



361 ] . Considering of the asymptotic behavior of J u (z) and Y u (z) around the origin and 



large \z\ we have the solution of Eq. ()2U[) to be of the form 

t = CJ v (z)+DY v (z), (21) 

where C and D are constants. The radial eigenfunction for the positive energy is 

u{p) = z v (CJ v (z) + DY u (z)). (22) 

The related probability density function of the above eigenfunction does not possess of the property 
of localization near the origin. 
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B. Three-dimensional system 

We use the same form of mass function fi(r) = r a as that for the two-dimensional system, then 
the radial Schrodinger equation for the three-dimensional quasi-free particle is 

ld 2 u(r) 12 -adu(r) 1(1 + 1) , . 

where u(r) is the radial wave function which has the relation to R(r) in Eq.© 

u {r) = ^1. (24) 

For the state of I = 0, Ea.(|23j) reads 

^ + 2 - adu ^ = 2\EKu(r). (25) 
dr z r dr 

Making the replacements of the following forms 



z = him r ^)/^ „ = ^ u{p) = , m ( 26 ) 

a + 2 2 + a 

we will turn Ea.(|25|) into the standard modified Bessel equation as Ea.(|15[). Now, v is restricted 
in the range (—5, 1) due to a > 0. 

Similar to the analysis of the two-dimensional system and with the requirement that u is finite 
when r is large, the general solution of Ea. (|15|) is 

Z = cK v (z), (27) 

where c is an arbitrary constant. The normalized solution of Ea.l|25j) is 

3(1-2") 



-Ci-u) 



(2\E\)^ l ~^ z v K u (z) 



(28) 



8\E\ sm( vir)( 1 - v) - , ? , , /— — 3 
_LJ \ — A ^^^-^(-(l - ^) v / 2|^|r 2 ( 1 —)). 

3^77 3 

The probability of finding the particle in the range r to r + dr is 

W(r)dr= / d6 sin c/? dysr 2 dr\ip\ 2 
io JO 

2-7T /*7T 

d0 / sin(^(iv?r 2 (ir|n(r)| 2 |yoo(6 l ,^)| 2 = |u(r)| 2 r 2 dr (29) 
Jo 

8\E\ sin(i^7r)(l — v) 2^ „ 2. , — r ? x , 

=-!— ! 1 — -ri— Kl(-(l - u)y / 2\E\r^—))dr. 

For the three-dimensional system, the relationship between the probability density W(r) and r is 
similar to that for the two-dimensional system. There also exists localization of the wave function 
for the s state around the origin. 
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Similar to the discussion for the two-dimensional system in the last subsection, the eigenfunction 
and the related probability density function for the positive energy in three dimensions have also 
no localization close the origin. 

IV. CONCLUSIONS AND DISCUSSIONS 

In this article, we discussed the properties of the s-state of a quasi-free particle with PDM. 
When mass varies with radial coordinate r, it induces additional effective potential which may 
produce central attractive force and may cause the wave function of the s-state to localize around 
the origin. In D, except D = 1, dimensional system with the mass of the particle being position- 
dependent, the localization appears, which contrasts sharply with the system of constant mass. We 
solve the radial Schrodinger equation to get the exact s-wave function for a specific mass function in 
two and three dimensions, respectively. Their corresponding probability densities and the related 
asymptotic behaviors are also analyzed, and the results are consistent with the general discussions. 

There are many kinds of mass functions that are used in the investigation of the systems with 
PDM, but it is difficult to get the exact solutions of the radial Schrodinger equations for these 
mass functions in higher dimensional systems. So some numerical studies are required for the 
understanding of the properties of the systems with these mass functions. 
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